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We analyze the relationship between the lowest energy configurations of an infinite harmonic 
chain of particles in a periodic potential and the evolution of characteristics in a periodically-forced 
inviscid Burgers equation. The shock discontinuities in the the Burgers evolution arise from ther- 
modynamical considerations and play an important role as they separate out flows related to lowest 
energy configurations from those associated with higher energies. We study in detail the exactly 
solvable case of an external potential consisting of parabolic segments, and calculate analytically 
the lowest energy configurations, as well as excited states containing discommensurations. 
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I. INTRODUCTION 

The Frenkel-Kontorova (FK) model is a classical infi- 
nite chain of atoms linked by elastic springs with equi- 
librium spacing /x, subject to an external periodic po- 
tential of period 2a^. It is a simple model for the de- 
scription of dislocations in solids but has also been ap- 
plied to interfacial slip in earthquakes^"—, the dynamics 
of DNA denaturation^, as well as other areasS. The FK 
model is characterized by the competition of two differ- 
ent length scales, /i and 2a, and two energy scales set 
by the external potential and the spring elastic energy. 
The static configurations of lowest energy have a rather 
complex dependence on these parameters: they can be 
commensurate or incommensurate with the period of the 
external potential^ and the transition between them is 
a critical phenomenon exhibiting scaling^^^^. The FK 
model has been analyzed based on a connection with two 
dimensional area-preserving maps, such as the standard 
mapi. Aubryi^ and Mather— have shown that the low- 
est energy configurations correspond to invariant sets of 
the associated maps that can be KAM tori. Cantor sets 
(Cantori), or limit cycles corresponding respectively to 
unpinned incommensurate, pinned incommensurate and 
pinned commensurate configurations. 

It was recently discovered that a continuum hydrody- 
namic type of evolution underlies the FK models^ii^. 
For the case of an elastic chain of particles embedded 
in an external potential, this evolution is governed by 
a periodically-forced inviscid Burgers equation and the 
associated flow of characteristics turns out to be closely 
related to the particle configurations. This connection 
was further developed by E and Sobolevskiii^iii (see also 
the review by Bee and Khanin~). 

The purpose of the present article is to illustrate and 
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further investigate the relation between FK models and 
its description in terms of a forced Burgers equation by 
explicitly working out an example. This is desirable 
for two reasons: On the one hand, the results obtained 
iniir— are mathematical, centering mostly around exis- 
tence theorems and properties of the flow with less em- 
phasis on connections with particular physical models. It 
would therefore be useful to consider an example that can 
be exactly solved using this approach, and thereby illus- 
trate explicitly how the flow properties relate to physical 
properties of the FK model such as the lowest energy con- 
figurations and excited states. On the other hand, most 
of the results on FK models focus on the lowest energy 
configuration from which other properties such as the 
stability region of a given configuration can be obtained. 
However, the Burgers description and in particular the 
flow patterns provide additional insight, allowing us to 
construct analytically the excited states with and with- 
out discommensurations and to see how the flow pattern 
changes as the parameters of the model are varied, par- 
ticularly near the boundary of a region of stability. 

The description of FK models by an underlying one 
parameter continuum fiow also constitutes a novel tech- 
nique which should be applicable to the larger class of 
problems^Tii mentioned above, as well as FK models with 
more complex external potential o^^'^° . The article is or- 
ganized as follows. In Section |IT] we briefly review the 
basic results due to Aubry, Aubry and collaborators, and 
Matherii^ii^i^. We then proceed in Section Hill to give 
a simple derivation of the relation between FK models 
and the periodically forced Burgers equation and discuss 
general properties of the solutions. In Section IIVI we 
focus on the FK model with a piece-wise parabolic po- 
tential. This model is interesting in its own right, as 
it is related to other FK models in the limit of strong 
external potential where the particles are conflned to re- 
gions close to the potential minima. We conclude with 
a discussion of our results and possible generalizations. 
Necessary background material as well as details of some 
of the calculations have been gathered in the appendix. 
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II. INFINITE HARMONIC CHAINS IN 
PERIODIC POTENTIALS 

We consider the static lowest energy configurations of 
an infinite cliain of particles of unit mass connected by 
Hookean springs with equilibrium spacing /i and subject 
to a periodic external potential with period 2a. The en- 
ergy of a configuration {yi} is given by the Hamiltonian 



Vi- li) + V{yi) 



, (1) 



is a Devil's staircase. Aubry has also shown that the 
lowest energy configurations {y^} are of the form 



Vi — /(*^ + a) = i£ + a + g{ii + a), 



(6) 



where g is periodic with period 2a and the choice of a 
only serves to determine which particle on the infinite 
chain is to be denoted the zeroth particle. The function 
g, or equivalently /, is called the hull function^. In the 
presence of KAM tori the hull function is continuous, 
whereas for a given i and sufficiently large r it becomes 
discontinuous. 



where for reasons that will become apparent soon, we 
have written the spring-constant as l/r. 

Any equilibrium configuration must satisfy the set of 
coupled difference equations 



Vi+i - 2yi + Vi-i - TV'{yi) = 0, 



(2) 



but solutions to the above equation will in general not be 
lowest energy configurations. 

Notice that due to the periodicity of the external po- 
tential V , if {yi\ is an equilibrium configuration, so is 
the configuration {yi -\- 2a}, where each particle has been 
shifted by an amount of 2a. Making the change of vari- 
ables y^ = y^ mod 2a and Pi+i = yi+i - y^, Eq. © 
becomes 



Vi+i = Ui 
Pi+i = Pi 



-Pi+i, 



(3) 
(4) 



which is a 2d Hamiltonian mapping on the cylinder S xTZ 
called a twist map. In the case of a sinusoidal external 
potential V{y) = 1 — cosijry/a), this map is also known 
as the standard map^^. 

Aubry and Mather have independently shown that the 
lowest energy configurations correspond to special invari- 
ant sets of the twist map, Eqs. with limit cycles 
corresponding to commensurate structures, whereas tra- 
jectories that are dense on KAM tori correspond to in- 
commensurate structures^. The KAM theorem applied 
to twist maps indicates that there is a critical threshold 
Tc such that for r > Tc all KAM tori have broken up 
and the incommensurate structures are dense on Can- 
tor sets (Cantori), that have measure zero. Thus for 
T > Tc almost all lowest energy structures are commen- 
surate (strong pinning limit). 



Aubry has shown that under general conditions 



to 



each lowest energy configuration {yi\ there is associated 
a winding number £ given by 



lini 



Vi+N — Vi+N' 

N-N' 



(5) 



which is the average distance between two neighboring 
particles in the lowest energy configuration. Moreover, £ 
as a function of the parameters /i and r takes constant 
values for each rational value of £/2a. Hence this function 



III. FK MODELS AND CHARACTERISTIC 
FLOWS IN A PERIODICALLY-FORCED 
INVISCID BURGERS EQUATION 

The determination of the lowest energy configurations 
as sketched in the previous section is rather indirect. The 
equilibrium equations Eq. or Eqs. ©-(HI), do not di- 
rectly yield them. Rather, the lowest energy configura- 
tions turn out to correspond to a highly special subset of 
orbits of the mapping. Their properties can be described 
qualitatively from the general properties of twist maps. 
Analytical calculations for a particular model are hard, 
and numerical simulations are difficult, since the trajec- 
tories in question are hyperbolic and hence numerically 
unstable^i^i. It would therefore be desirable to obtain 
such configurations in a more direct way analytically. 



A. The Fundamental Catastrophe 

Consider a particle in an arbitrary potential (not nec- 
essarily periodic), as shown in Fig. [1] to which there is 
connected a spring of spring constant l/r. Let xi de- 
note the location of the particle, while X2 denotes the 
position of the endpoint of the spring, such that xi = X2 
corresponds to the spring being unstretched. 

We are interested in the position of the particle xi as 
a function of X2 . This can be found from minimizing the 
Hamiltonian 



H{x2,Xi) 



1 

27 



{X2-Xif +V{xi) 



with respect to xi so that 



Xi = argimn H {x2, Xi) 



(7) 



(8) 



The problem is non-trivial due to the non-convexity of 
V{x). Differentiating H{x2,xi) with respect to xi, 



X2 — xi + tV'{xi). 



(9) 



For T sufficiently large this equation does not necessarily 
have a unique solution for all values of X2 anymore, as 
the dotted curve in Fig. [^fa) shows. In order to obtain 
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FIG. 1: Particle in an external potential V{x) connected to 
a spring with stifltness 1/r. The positions of the particle and 
the endpoint of the spring are given by xi and X2, respec- 
tively, such that xi = X2 corresponds to the spring being 
unstretched. 



a single- valued dependence of xi on X2 an additional as- 
sumption is needed: we require that the woA done in 
moving the end point of the spring is equal to the change 
in total internal energy of the particle. The latter is given 

by 



Hint{x2,T) = H{x2,Xi)\^^-_ 



-■Xi{x2,t) 



(10) 



resulting in the red curve in Fig.[2lja) with a jump discon- 
tinuity which corresponds to the particle abruptly switch- 
ing wells as X2 is increased. As a result, the jump of the 
particle from one well to the other does not generate any 
heat and the process is adiabatic. The assumption made 
is thus thermodynamical. The work done on the system 
is readily shown to be continuous in X2, and Hint{x2,T) 
is continuous in X2 as well. In terms of the location of the 
jump discontinuity this implies that the areas bounded 
by the dashed curve and the discontinuity in Fig. [IJb) , 
have to be equal. This is the familiar Maxwell equal-area 
construction. 

Let us now reconsider the internal energy of the system 
by treating the reciprocal spring constant r as an addi- 
tional variable. By definition, Hi^t {^2 ,t) is a potential 
so that the associated force 



Fix2,T) 



dH-n 
dx2 



(11) 



corresponding to the restoring force at the end point of 
the spring, must be conservative. From simple mechani- 
cal considerations it is also clear that 



F{x2,t) 



dV 
dxi 



(12) 



where Xi satisfies Eq. ([9)). Note in particular, that for 
T = 0, corresponding to an infinitely stiff spring, we have 
xi = X2, SO that 



dV 
dxi 

Combining the above equations, we find. 



(13) 



F{x2,t) = F{xi,0), where X2 = xi -tF{xi,0). (U) 

In other words, F{x, t) remains constant on the line x = 
xi — tF[xx,Q) with < i < T. Physically, this is just a 
restatement of the fact that the forces on both ends of 
the spring are equal. Mathematically, however Eq. (fT4|) 
implies that F{x, t) is a solution of the inviscid Burgers 
equation 



dt ox 



(15) 



with initial condition F{x^ 0) = —dV/dx. 

It is more convenient to write the inviscid Burgers 
equation in its more familiar form by letting u{x, t) = 
—F{x,t). Denoting partial derivatives by subscripts, we 
have 

Ut + uux = (16) 
with u{x,Q) = dV/dx and from Eq. ([M]) . 

u{x, t) u{x(i, 0), (17) 
for x, xq and r satisfying the characteristic equation 



X = xo + Tu{xo, 0), 



(18) 



which are lines on the xt plane. 

Although u is constant on the characteristics, whenever 
du/dx < 0, characteristic lines intersect, corresponding 
to multiple- valued solutions. As we have shown, this situ- 
ation is resolved by the equal-area construction and gives 
rise to a discontinuity, a shock. Such solutions are called 
weak since, Eq. (jl6p can only be satisfied in a weak sense, 
due to the discontinuities of u. Further relevant details 
on weak solutions of the inviscid Burgers are provided in 
Appendix [X] 

It is important to stress that the weak solutions are not 
a mathematical artifact, but follow from thermodynam- 
ical considerations. To see this more clearly, consider 
the particle-spring system embedded in a heat bath of 
temperature T. The position of the particle xi can be 
thought of as an internal variable and we consider the 
partition function 



\2ttt 



1/2 



e-^i—r e-^^(^') dx' 



(19) 

as a function of the external variable x, with e(x, r) be- 
ing the corresponding free energy. An x-independent 
pre-factor can be arbitrarily chosen, but with the choice 



4 



(a) 



(b) 




(c) 




FIG. 2: (Color online) (a) The Lagrange map yielding the position xi of the particle in Fig. [T] as a function of the endpoint 
of the spring X2- At a certain X2 the particle abruptly jumps wells. The dotted line is a plot of Eq. The actual location 
of the discontinuity on the red curve is obtained from an adiabatic condition, (b) The derivative of the internal energy of the 
particle Hint as a function of the endpoint X2 of the string. The dotted curve is the multiple-valued solution Hl^t containing 
overhangs. The location of the discontinuity X2 is such that the areas bounded by the dotted curves and are equal, (c) The 
internal energy Hint along the corresponding two paths in (a). 



made above e"^"^*^^'*^ is a solution of the diffusion equa- 
tion with diffusion constant kT/2. The Cole-Hopf 
transformatio n^^'^^ u(x,t) = ex{x,t) yields the viscid 
Burgers Equation for u: 



Ut + UUx 



kT 



(20) 



with initial condition w(a;,0) = V'{x). The same weak 
solution also follow from the solution of Eq. (|20l) at non- 
zero T in the limit T — >■ 021"—. Thermodynamically, this 
corresponds to cooling the system quasi-statically to zero 
temperature resulting in a lowest energy configuration, 
c/. Fig. [2][c). The internal energy Hi^t{x,T) is thus the 
zero-temperature limit of the free energy e(a;,T). 



B. Burgers Evolution of the FK Model 

The construction presented in the previous section can 
be utilized to treat the lowest energy configurations of 
the FK model with Hamiltonian Eq. ([T|). Consider a 
semi- infinite chain with particle configurations {yi\ such 
that —CG < i < n. Denote by H^"^(yn,T) the energy of 
a lowest energy configuration of the semi-infinite chain 
with its end point fixed at i/n- Owing to the periodicity 
of the external potential V{y) = V{y + 2a), Hi^^{y„,T) 
must also be periodic in i/„ with the same period 2a. 
Now add another particle to the right end of the chain. 
The total internal energy H^"^^^^ (?/„+i, r) of the resulting 

in). 



chain is related to H^^^{yn, t) as 



H.. 



[n+l] 



{yn+l,T) = V{yn+l) 



■ mm ' 



n+l 



Vn - M) 



2t 



(21) 



Recalling that /i is the equilibrium spacing of the 
springs, we can make a change of coordinates to posi- 
tions relative to the endpoint of each unstretched spring 
as 



Xi = yi - iM, 
and Eq. becomes^^i^ 

i7il:;+^^(x„+l,T) = V{Xn+l + (n + 

I • I (^n+1 ^ Xn) Tj(n) / \ 

™ I 27 



(22) 



(23) 



In this form the above equation closely resembles the 
problem presented in the previous section, Eqs. ([7]) and 
(|5]), and we can evaluate the expression to be minimized 
on the RHS of Eq. via evolution of the inviscid Burg- 
ers equation, as follows: 

Let 



Un{x) 



dx 



(24) 



then following the steps of the previous section and treat- 
ing f as a variable, we next define 



H-Ux.t) = min |1 {x - x'f + i?iL"^x',T) 



(25) 



Hence u(a;, t) = dHint [x, t) /dx satisfies the inviscid Burg- 
ers equation 

ut + uu.:^ = Q for < i < T (26) 
with u{x, 0) = Un{x), so that from Eqs. ([23]) and (HJ]) 

Un+i{x) = V {x + {n + + u{x,t). (27) 
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Moreover, the position x„ of the n particle as a func- 
tion of the position Xn+i of particle n + 1 is given by the 
characteristic mapping 

Xn+l ^ Xn+ Un{Xn)T. (28) 

The relations Eqs. (f24|) - (|27|) actually prescribe the 
evolution of u under a forced Burgers equation. Defining 

"(a;,Olt=(„r)+ = ""(2;)' (29) 
the evolution equations become 

Ut + uUx = 0, nr < t < (n + l)r (30) 
w(a;,Olt=(„r)+ = "(a;,i)lt=(„r)- + V'{x + n^i), (31) 

which is equivalent to the periodically forced Burgers 
equation: 

oo 

Ut + UUx = ^ 5(t — tit) V'{x + nfi), (32) 

n=0 

with initial condition u{x, 0^) — 0'^^'. 

The flow of characteristics, Eq. (1^51) . under forced 
Burgers evolution implicitly defines the characteristic 
backwards map (also known as the Lagrange map) , such 
that, given a final time tg, for all t < to 

x{t) = x{to;t). (33) 

For the configurations {xi} of the semi- infinite chain with 
the outmost particle n being at a;„, this implies that for 
all i < n 

Xi = x{nT;iT), with Xn=x{nT;nT). (34) 

We have therefore shown that a continuous one- 
parameter flow embodied by the Lagrange map Eq. p3p , 
underlies the equilibrium configurations Eq. p4p of the 
discrete mass-spring system. The Lagrange map in turn 
is given by the backwards flow of the characteristic tra- 
jectories of the forced Burgers evolution Eq. ([32]). Within 
this description the time-like evolution parameter i is a 
material coordinate corresponding to the building up of 
springs by the continuous addition of material with elas- 
tic modulus /i/r. 

Our derivation of the connection between the discrete 
particle configurations of a harmonic chain of particles 
and the characteristic fiow of a forced Burgers equation 
has been based on an analysis of the forces acting on the 
particles. The connection between a general class of dis- 
crete minimization problems such as Eq. (|2ip and certain 
one parameter flows was established first independently 
by Jausslin, Kreiss and Moser-i^ and E, Khanin, Mazel 
and Sinaii^, using variational methods. The connection 
with FK models in particular was developed further by 
E and Sobolevskiii^. 



C. Properties of the Characteristic Flow Pattern 

For the FK model the relevant results are as 
foUowsi^"— : (i) To each asymptotic solution u{x,t) 
there corresponds a fiow pattern of characteristics 7 
with 7(^0) = which are traced backwards in time, 
t G (-00, to], cf. Eqs. (1321) and These charac- 

teristics 7 cannot cross each other and by construction, 
will never terminate in a shock. They are called one- 
sided minimizers. We should re-emphasize that by def- 
inition one-sided minimizers flow backwards in time. In 
the case of FK models they generate the lowest energy 
particle configurations of a semi-infinite chain with the 
outermost particle fixed at xq. (ii) Among the one sided 
minimizers there exists a subset of minimizers that have 
the additional property that when traced forward in time, 
t > to, they never merge with a shock. These minimiz- 
ers are the global minimizers and they correspond to the 
lowest energy configurations of the bi-infinite chain. As 
t — > ~oo, the one-sided minimizers converge to one of 
the global minimizers. (iii) Given a time t, the set of 
points Xs such that minimizers immediately to its right 
and left {xf and xj ) converge to different global minimiz- 
ers, constitute the locations of global shocks. Thus global 
shocks, if present, have the property that they can never 
disappear, as they separate the fiows of one-sided min- 
imizers that approach different global minimizers. (iv) 
All minimizers associated with an asymptotic solution 
u{x,t) have the same winding number i/2a, correspond- 
ing to the average spacing of particles of a configuration, 
cf. Eq. ([5]). (v) Pinned particle configurations are char- 
acterized by the presence of shocks in the flow patterns. 
For rational £/2a the flow pattern turns out to always 
contain shocks and the particle configurations will thus 
be pinned. For irrational £/2a, depending on the external 
potential, the asymptotic flow pattern may or may not 
contain any shocks. These cases correspond to pinned 
and sliding incommensurate configurations, respectively. 

For the particular FK model with piece- wise parabolic 
potential, which is the case we will be concerned here, 
the external forcing always contains a shock. Since a 
shock once present cannot disappear but at most will 
merge with another shock, shocks will be present in the 
fiow pattern and the resulting particle configurations are 
pinned. Furthermore, almost all conflgurations (except a 
subset of measure zero) turn out to have rational winding 
numbers i/2a = r/s. The fiow pattern is periodic in time 
t with period sr— and thus global minimizers corre- 
spond to a periodic configuration of particles with period 
s. It is not hard to see that there must be s global mini- 
mizers: At any time t = nr their locations x correspond 
to the s distinct locations of particles in the periodic low- 
est energy configuration. Since starting from a time ^Qj 
the backwards flow of one-sided minimizers converges to 
a global minimizer, these must converge to one of the 
s global minimizers. Given that the configuration space 
X is periodic (with the period of the external potential 
2a), the unit cell must contain s global shocks separat- 



V'{x)k 



FIG. 4: The profile u{x, t) and the parameterization of its 
linear segments. 











y2a 




/2a *'x 



FIG. 3: Top: The picccwisc parabolic potential V{x), 
Eq. (|35[) . Bottom: The corresponding piecewise linear pro- 
file V'{x). 



ing the backwards flows of one-sided minimizers towards 
their associated global minimizer. In particular, the flow 
of the minimizers in the xt plane will be confined to the 
interior of s strips that are bounded by the trajectories 
of the global shocks and that each contain a global min- 
imizer. The (backwards) flow of one-sided minimizers 
remains thus inside their respective strips and thereby 
converges towards the associated global minimizer. 



IV. BURGERS DESCRIPTION OF THE FK 
MODEL WITH PARABOLIC POTENTIAL 

In this section we calculate the flow patterns of a FK 
model with a piece-wise parabolic potential. This case 
also corresponds to a strong pinning limit in which the 
external potential is so strong that particles in the low- 
est energy configuration are confined to the vicinity of 
the minima of the potential wells, that can be treated 
approximately as parabolic. The lowest energy configu- 
rations of this chain were calculated exactly by Aubry^^. 
The purpose of this section is to recover Aubry's solutions 
and to demonstrate how the forced Burgers evolution ap- 
proach yields additional results and insights. 



A. Parameterization of the Burgers Profile and its 
Evolution 



The external potential is 



V{x) = -Ao 



X — 2a Int 



X + a 
2a 



(35) 



where 2a and Aq > are the period and the strength 
of the potential, respectively. We consider a unit cell 
that extends from —a to a. The forcing of the Burgers 
equation, Eq. p2p . is given by V'{x) which is a series 
of ramps as shown in the bottom part of Fig. [3] The 
x-intercepts correspond to the potential minima. The 
continuity of V{x) across the boundaries of the unit cell 
further implies that the total area under the profile from 
—a to a is zero (area constraint). 

It is not difficult to see (c/. Appendix]^ that evolution 
under Eq. ([32]) is such that for any time t, the Burgers 
profile within a unit cell consists of parallel straight line 
segments of slope \{t) terminated by shocks. A sample 
profile is shown in Fig. SI In what follows we will be 
making use of certain facts about the evolution of such 
profiles. The relevant results have been derived in Ap- 
pendix lA 11 

At any instant of its evolution, the profile of u{x, t) is 
completely determined by a set of parameters^: Each 
line segment is part of an infinite line of slope \{t). Since 
the segments are confined between shocks, the position 
of the shocks determine the intervals that the segments 
occupy on their respective lines. Along with their slope, 
these lines are determined by their x-intercepts. Hence, 
if there are k shocks inside the unit cell, there are k -I- 1 
segments. Numbering the segments within a unit cell 
from left to right as to k, we will denote the right ter- 
minations of each segment as and the corresponding 
x-intercept by v^^\ Note that due to the periodicity 
across the unit cell, the intercept associated with the last 
segment is given as i^*-**^ — 2a + v'^'^\ Including the slope, 
2k + 1 parameters are required to determine the profile 
u{x,t) completely. 

The evolution of u{x, t) under Eq. p2|) consists of two 
parts: the evolution step Eq. ([50)1 when a new spring is 
added to the end of the chain, and the particle insertion 
step Eq. (PT|) . where a new shock is inserted, c/. Fig. ([SJ. 
During the evolution step the slopes of the segments flat- 
ten according to Eq. (jA12l) and the shocks move and 
merge upon collision. Whenever a new shock is inserted, 
one of the linear segments of the profile will be split into 
two by the shock discontinuity of V'{x + n^) (unless it 
happens to coincide with the boundary of a segment), 
and the slopes of the profile will be incremented by Aq, 
cf. Eq. (1211). 
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FIG. 5: Top: The profile m and the profile V' to be superposed 
(dashed). Bottom: the resulting profile after the addition of 
the two profiles. 



Let us first consider the evolution of X{t). Denoting 
by the profile slope just before and after addition of 
a particle at time nr, we see from Eq. (|A12p that 



A, 



An 



A^ 



1 + XT^T 

This recursion converges to a stable fixed point 



A* 



1 



tAo 1 



1 



4 
tXq 



(36) 



(37) 



In what follows we will also need Al/A!j_, which, noting 
that Al = Al — Aq, turns out to be 



X*_ 

Al 



1 + A* 



tAq 
2 



tAo 
2 



4 
tAo 



77, (38) 



coinciding with 77 irt^i. Asymptotically, the profile slopes 
right before and after a particle addition are thus given 
by A j_ . We will henceforth assume that sufficiently many 
particles have been added to the chain that the profile 
has reached its asymptotic slope. 

During the evolution step a segment k disappears 
whenever shocks ^(''""i) and ^'^'"'^ merge. The evolu- 
tion of a segment k that survives the evolution step 
nr <t < {n+ 1)t is given as, Eq. (|A11[) . 



{x,t) 



X* 



1 + (t - nT)Xl 



(39) 



for ^('=-i)(i) <x < ^'•''Ht). Note that z/C') is constant for 
segments k that survive the evolution step. The location 
of the intercept points v^''^ will generally change during 
the particle insertion step, since besides creating a new 
segment, the slopes of all segments are augmented by 
Aq, while the locations of the shocks already present 
remain unchanged. 



The mapping of j^^'^^ during particle insertion can be 
worked out and is illustrated in Fig. [5] Denote by 
the location of the segment intercepts before and after 
the particle addition, and let b be the location of the 
intercept of the left segment of the new shock to be added 
(see Fig. [5]). The new shock thus is at ^('1'=*) = b + a and 
the intercept to its right is at & + 2a. One finds that 



where 77 is defined in Eq. (|38|) and 



2a, 



^(fc) ^ 



(40) 



(41) 



Note that the segment intercepts, Eq. (PH)) . are attracted 
toward their respective b^'^^'s, since rj < 1 unless Aq = 0. 

If, as will generally be the case, the new shock splits 
a particular segment k into two, this will also create an 
additional intercept, that is given as 



(new) 



■qiy--' + (1 - 77)6. (42) 
The ordered list of intercepts after insertion is thus 

(43) 



7.(0)+, 



(fc-l)+ („cw)+ + 



Since 7y('^)+ remains constant for segments surviving the 
evolution step, we will henceforth drop the + super- 
scripts. 

We consider next the evolution of the shocks. De- 
note by S,^^^ and i^'-'^-' the location of the shock and the 
corresponding segment intercept, respectively, right af- 
ter a shock insertion. As we show in Appendix lA II the 
shocks move at constant velocity w^*^) and we find, using 
Eq. 1MB, 



,{k) 



,{k+l) 



(44) 



The final location Q of a shock that survives the evo- 
lution step without merging with another shock is then 
found using Eq. ([55]) . as 



(fc) 



1 



(fe) 



1-77 T/^'') + 



(45) 



Finally, the parameter h indicating the location of the 
zero intercept of V'{x-\-{n-\-l)^) in the co- moving frame, 
evolves according to 



(46) 



Let subscripts j denote the times t — (jr)^ right after 
shock insertion. Together with the rules of how to handle 
colliding shocks given in Appendix lA H we thus have a 



discrete dynamical system for the variables v^'^'^ and ^('^' 



that underlies the evolution of u{x,t), Eqs. (|B1|) - (|B5|) 

For the FK model with piece- wise parabolic potentials, 
all lowest energy configurations are commensurate with 



an average spacing Ijla = r/s, where r and s are rel- 
atively prime integers. Correspondingly, the asymptotic 
behavior of u{x, t) is periodic up to a shift in the sense 
that, for all x and t 



u^{x,t + st) = u^{x + Sfl,t). 



(47) 



Focusing on the moment right after a particle addition, 
this means in particular that there are s different profiles, 
w+(x,0),w+(a;,r),u+(a;, 2r), . . .,u+{x, (s-1)t) that turn 
out to be related to the semi-infinite chain with its end 
particle located at one of the topologically distinct s lo- 
cations of the lowest energy configuration, as we will see 
shortly. 

Before proceeding with an analytical derivation of the 
steady-state profiles and the associated characteristic 
flow patterns, it is instructive to look at some steady- 
state solutions obtained by numerically evolving the pro- 
file parameters j/ and ^. 



B. Steady-State Profiles and Flow Patterns - 
Overview 

Figure |6] shows the shock trajectories associated with 
the steady-state flow of the forced Burgers equation with 
parameters ///2a = 0.8495 and Aq = 0.4, corresponding 
to i/2a = 1. Since we will be only interested in steady- 
state flow patterns and not the transients, we have reset 
time to t — 0. From Eq. (|47|) we see that the profiles 
u{x,t) at i = nr and t = (n + 1)t are equivalent up to 
shifts by —fj,. At each time t — nr a new shock is inserted 
due to the particle addition step and the insertion point 
is marked by a triangle. The bottom figure shows u{x, t) 
at t = 0+. 

Recall that with respect to the co-moving frame the 
unit cell is moving by an amount of — /i from one parti- 
cle addition to the next, c/. Eq. ([22l) . The boundaries 
of the unit cell coincide with the shock insertion loca- 
tions, corresponding to the cusps of the external potential 
Eq. p5|) . They therefore also mark the location of the 
unit cell y G [—a, a) in the co-moving frame, correspond- 
ing to y — ±a in the unit cell coordinates. The locations 
of the minima of the external potential lie half-way be- 
tween these two points at y = in the cell reference frame 
and are marked by open circles. In terms of coordinates 
y we have strict periodicity, u^(y,t + r) = UH.(y,t), i.e. 
without shift. 

Note how a newly inserted shock can survive subse- 
quent insertions thereby creating a tree-like structure. 
We will refer to such structures as shock-trees. In fact, 
each inserted shock eventually collides with a previously 
inserted shock and the corresponding collision events 
have been indicated in the figure by red asterisks. Due to 
the periodicity of the shock configurations, the life-time 
5tc of a newly inserted shock is constant. In the figure 
shown, dtc is between 5t and 6t implying that at any 
insertion time t = {nr)^ there are 7 shocks including the 
newly inserted shock (two of the shocks are too close to 




FIG. 6: (Color online) Top: Shock trajectories of the steady- 
state flow for parameter values fj./2a = 0.8495 and Aq = 0.4. 
The period in this case is r, corresponding to £/2a = 1, and 
thus the profiles of u(x,t) at t — nr and t — {n + l)r are 
equivalent up to shifts by —fi, cf. Eq. (|47p . At each time 
t = nr a new shock is inserted due to the addition step and 
the insertion point is marked by a triangle. The red open 
circles show the locations of the well minima of the external 
potential (see text for further details). Bottom: Profile of 
u{x,t) at t = O""". 



be discerned, see bottom half of Fig. [6]) . It is not hard 
to see that the time-periodicity of the profile also implies 
that all branches of the shock tree starting with a newly 
inserted shock are identical. They are merely nested on 
the tree due to their different creation times. 

Once the inserted shock merges with a previously in- 
serted one, this shock proceeds to evolve until the next 
collision. One can think of the latter as an ever-present 
shock that keeps on absorbing newly inserted shocks, 
thereby constituting the trunk of the shock tree. This 
is the global shock introduced in section IIII C[ while 
the shocks constituting the branches are the secondary 
shocks^^ . 

Figure [7] shows the shock trajectories of the steady 
state flow for £/2a — 2/5. There are 5 shock trees. For 
each tree there is an insertion time at which a new shock 
is added. By the periodicity of the flow it follows that 
shocks are added to a given tree periodically (the period 
in this case being 5t) . Observe the "feeding order" of the 
trees. The immediate right neighbor of a tree on which a 
shock has been inserted is "fed" at the second subsequent 
insertion. As has been shown^^'^^ {see Section llll CI) , the 
shock pattern for a conflguration with £/2a — r/s, with 
r and s irreducible integers, will always contain s shock 
trees and the time periodicity of the pattern will be t = 
ST. Labeling the trees from left to right as 0, 1, 2, . . . , s — 
I, the feeding order of the trees turns out to be given by 
subtractions of r mod s, as is proven in Appendix IB II 
From the periodicity sr of the flow pattern it also follows 
that for a given shock tree all its secondary shocks are 
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FIG. 7: (Color online) Shock trajectories for fj,/2a = 0.39 
and Ao = 0.2 at steady-state. The equilibrium spacing of 
the particles is I /2a = 2/5. Note that the figure contains 
s = 5 shock trees and that the tree patterns are periodic 
in t = 5r up to an overall shift. For any shock tree, the 
secondary shocks (light blue) are either all to the left or the 
right of their global shock (dark blue). It turns out that at 
their respective insertion times the trees (from left to right) 
contain 2, 8, 7, 8 and 2 shocks, but most of these are too close 
to the global shock to be discerned. 



FIG. 8: Shock trajectories and backwards flow of the char- 
acteristics (minimizers) for /i/2a = 0.8495 and Ao = 0.4 
at steady-state. The average spacing of the particles is 
I /2a = 1/1. The annotation of the shocks and unit cell 
boundaries are as in Fig. [G] The trajectories in green denote 
one-sided minimizers, whereas the dark green trajectory is the 
global minimizer that generates the lowest energy configura- 
tion in the co-moving frame. The locations corresponding to 
the positions of the particles in this configuration have been 
marked on the global minimizer by black solid circles. 



either to the left or right of its global shock. We will 
refer to such trees as left and right trees, respectively. In 
other words, the time periodicity of sr of the flow pattern 
and the presence of s shock trees onto each of which a 
single secondary shock is inserted during the period st, 
implies that a shock tree cannot have branches on both 
sides of its global shock. However, both types of trees 
can coexist, as seen in Figs. [7] and HI 

We now turn to the flow pattern associated with shock 
trajectories at steady state. The lowest energy configura- 
tions of a semi-infinite chain with its end point fixed are 
generated by the characteristic trajectories traced back- 
wards in time according to Eq. p4l) . These trajectories 
are the one-sided minimizer^. Since the flow at steady 
state is time-periodic with period st, it is sufficient to 
know the backwards flow for times t £ [0, st) and for all 
X € [—a, a), generating a map that can be iterated. 

We expect that as the characteristics are traced back- 
wards in time, corresponding to particles deeper and 
deeper inside the chain, the effect of the location of the 
particle at its end point will diminish. As we will show, 
the effect of the boundary decays as 77* . This means that 
as the characteristics are traced further back, the charac- 
teristic trajectory will approach a limiting cycle (due to 
the periodicity of the steady-state flow pattern). Since 
the effect of the boundary will have vanished, the par- 
ticle configuration generated by the limit-cycle must be 
the lowest energy configuration of the bi-infinite chain. 



This confirms the remarks in Section fill C[ namely that 
one-sided minimizers when traced backwards in time will 
converge to a global minimizer. Equivalently, the global 
minimizers when traced forward and backwards in time, 
generate the lowest energy configuration of the bi-infinite 
chain. 

Figure [8] shows the flow pattern associated with Fig.|6l 
The one-sided minimizers correspond to green lines, while 
the global minimizer is indicated by a darker green line. 
We have £/2a — 1/1, meaning that each potential well 
contains one particle. The location of the particles of the 
lowest energy configuration in the co-moving frame are 
shown as black solid circles and they necessarily lie on 
the global minimizer. Furthermore, they turn out to co- 
incide with the locations of the minima of the potential 
wells (red open circles) . Figure |9] shows the fiow pattern 
for £/2a = 1/7. Note how the one-sided minimizers when 
traced back in time flow onto one of the seven global min- 
imizers. Again, the region between any two shock trees 
contains precisely one global minimizer and the corre- 
sponding particle configurations in the fixed frame are 
all equivalent up to an overall cyclic permutation. 



C. The Fundamental Shock Tree 

Given a configuration with £/2a = r/s, the correspond- 
ing flow pattern will contain s shock trees and also s 
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global shocks. The global shocks have the property that 
they do not disappear. The secondary shocks constitute 
the branches of the shock tree that will eventually merge 
with the global shocks. The presence of a gap region in 
between shock trees that is bounded on each side by a 
shock, implies that the corresponding segment of u{x,t) 
will never disappear, since the shocks at its boundaries 
will never merge with each other. We will refer to these 
segments as global segments. Likewise, the intercepts as- 
sociated with such segments will never disappear and we 
will refer to them as the global intercepts. 

One note of caution is in order. Recall that our conven- 
tion has been to identify the intercept of each continuous 
segment of u{x,t) with the same label as the shock that 
bounds it on the right. With this convention one must 
be careful since the global shock bounding a global seg- 
ment may fall to the left of the segment, and therefore 
the global segment and the global shock may not have the 
same label. This occurs for a left tree: By definition, the 
rightmost shock in a left tree is its global shock. The seg- 
ment bounded on the right by this shock is however not a 
global segment, since this segment and its intercept will 
disappear when the secondary shock bounding the seg- 
ment from its left merges into the global shock. It is not 
difficult to see that for a left tree the segment bounded 
on the left by the global shock is a global segment. Thus 
its label will be that of the shock lying immediately to 
the right of the global shock. In general this shock will 
belong to another tree^. In the case of a right tree such 
a situation does not occur. The segment associated with 
the global shock is global as well and therefore carries 
the same label. 

The segments of u{x,t) associated with the global in- 
tercepts are by construction the segments of u that span 
the gap region in between two neighboring shock trees. 
This region contains the global minimizer and moreover 
governs also the backwards flow of one-sided minimizers 
in its vicinity. Except for the case where a right tree is 
immediately to the right of a left tree, the global region 
will be bounded on at least one side by secondary shocks. 
Thus whenever a new secondary shock is inserted, the 
corresponding global segment will be intersected, spawn- 
ing off two intercepts: one that remains global and one 
that is associated with the secondary shock. We will refer 
to the latter as secondary or local intercept. An example 
is shown in the right panel of Fig. (TU] Here the flow 
pattern contains a single shock tree that is of right type. 
Thus the global region is bounded by the two sides of 
the tree, which is equivalent to being bounded by two 
right trees. The global intercepts are shown in red, while 
the secondary intercepts are in green. The bifurcation 
into two intercepts occurs aX t — ut but for clarity has 
been offset in time by a small amount. As can be seen, 
the segment associated with the global region is bisected 
by successive shock insertions. The corresponding global 
intercept spawns off local intercepts that evolve on their 
own. The local intercepts are terminated when the corre- 
sponding secondary shock merges with the global shock. 




FIG. 9: (Color online) Shock trajectories and backwards flow 
of the minimizers for /i/2a = 0.155 and Ao — 0.2 at steady- 
state, corresponding to £/2a = 1/7. The flow pattern is 
marked by the presence of s = 7 shock trees. The anno- 
tation is as in Fig. [S] Note that for some of the trees the 
main shocks (dark blue) are obscured by the global minimiz- 
ers (dark green), since they are very close to each other (see 
text for further details). 



The time periodicity of the steady state fiow, Eq. (|47l) , 
implies that the number of shocks and intercepts is con- 
served across a period. Since a shock insertion always 
adds a new shock and a corresponding intercept, while 
a shock collision removes a shock and an intercept, it 
follows that during the period st the number of shocks 
inserted must equal to the number of shocks that merged 
with the global shocks. 



D. The Global Intercepts 

The global intercepts v can now be obtained as follows. 
By definition, global intercepts do not disappear and 
their locations remain constant during Burgers evolution. 
However during shock insertions they are remapped ac- 
cording to Eqs. dHU, (|B2|) and (|B3l) . This mapping in 
turn depends on whether the shock associated with v is 
to the left or right of the newly inserted shock. Thus 
we first have to determine the sequence bj*"' of Eq. (jB3P 
governing the evolution of the global intercepts. This is 
related to the "feeding order" in which new shocks are 
inserted into the trees and the details are given in Ap- 
pendix IB II For the segment associated with the global 
intercept of a shock tree the result is. 



bj = —jfi + 2a Int (j + S)- 



(48) 



where henceforth we will omit the superscript (fc), when- 
ever the segment in question is global. The integers 
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5 = 0, 1, . . . , s — 1 are each associated with one of the 
s shock trees. They are related to the relative time lags 
of each shock tree from their next shock insertion. 

By definition, the global intercepts survive all evolu- 
tion steps. This turns Eq. ((B2|) together with Eq. p8)) 
into a recursion relation that can be solved. For our pur- 
poses it is more convenient to consider the intercepts in 
the unit cell coordinates defined as follows: 

^3 = ^j~b]^ (49) 

which defines the location of the intercept relative to the 
location of the minimum of the potential well, so that 
Vj e [—a, a). As will be shown shortly, without loss of 
generality we will consider the shock tree for which (5 = 0. 
Then, the above definition along with Eqs. (|B2p and pS)) 
yields the following recursion 



where Xj = {bj+i — bj + /i)/2a is given as 
=Int((j + l)-) -Int(j- 



(50) 



(51) 



Note that Xj is periodic, Xj+s = Xj: well as xo = and 
Xs-i = 1 (except for the case r = s = 1, for which xo — 
1). In the cell coordinates, due to the time periodicity 
of the profile we also have I'j+s = i^j- Upon solving the 
recursion we find 



1^0 = Ai 



V 



2a 



s-l 



1 — 1] 1 — 77** 



E^'"^». (52) 



i=0 



The remaining i>j, for j = l,2,...,s — 1 can be then 
found as 



1-77 



(53) 



2a 



I ' 1=0 i=0 ) 



From the definition of Xj j Eq. (I5l1) , and its relation to 
bj, Eq. it is evident that a non-zero S will induce a 
cyclic shift of x by an amount of S. We thus define 



(54) 



Xi X{i+5) mod 



so that for 6 > this is a left shift. Using Eq. (|52p to 
define a function i>[x] such that i>o = v[x], the periodicity 
of i>j implies 



(55) 



At any time t — nr, the global intercepts associated 
with each of the s shock trees only differ by the time of 
last insertion of a shock, as captured by the distinct val- 
ues of 5. Thus from Eq. ([55)) we see that the set of global 
intercepts at any given time must also coincide with the 
set {I'j} (up to a cyclic permutation of its elements). In 



other words, {i/j} not only corresponds to the periodic se- 
quence of global intercepts associated with a given shock 
tree during its time evolution, it also corresponds to the 
set of all global intercepts associated with the s shock 
trees at any given time t — nr. Thus with respect to 
their global intercepts all shock trees are alike, differing 
only in their respective shock insertion times. 

In fact, it can be shown that essentially the same holds 
true for all intercepts associated with the shocks as well 
as the secondary shocks themselves: If we label the sec- 
ondary shocks and their corresponding intercepts at a 
time t — jr on a given tree by a superscript k, then 
it turns out there is an infinite sequence ^j"', i/j, with 
j = 0, 1, . . . , s — 1, and fc = 0, 1, 2, . . ., such that for each 
shock tree a with Kq, -f 1 secondary shocks, the actual 
secondary shocks and intercepts are subsequences termi- 
nated at k = Ka- This reduces the problem of obtaining 
the steady state shock pattern to finding the s global 
shocks and their k values. The calculations are rather 
involved and will be carried out elsewhere^. Instead, 
here we will restrict ourselves to s = 1 which is a special 
case of the above and already contains most of the rele- 
vant features of the general case. On the other hand, the 
global minimizers and thus the lowest energy configura- 
tions can be determined from the global intercepts alone, 
which we have just found for all r and s. We will carry 
this out next. 



E. Lowest Energy Configurations 

We turn first to the evolution of characteristics inside 
the global segments. As mentioned before, the charac- 
teristics x{t) associated with the lowest energy configu- 
ration are the global minimizers and have the property 
that they are periodic up to a shift. 



x{t + st) ~ x{t) — Sfl. 



(56) 



Our goal will be therefore to write the flow equation of 
characteristics inside the global segment and then impose 
the periodicity condition. Since we are interested in a 
periodic solution this calculation can be done forward or 
backwards in time, but it turns out to be more convenient 
to consider the forward flow of characteristics. Letting 
Xj = x{jT), the characteristic equation for any Xj within 
the global segment is given by 



Xj + l = Xj +X*^TiXj - Vj) 

_ 1 1 -?7 
7] rj ■'' 



(57) 



Introduce the unit cell coordinates i/j as 

Xj = y, + bj (58) 

so that, using Eqs and ([5T|) . the recursion for y be- 
comes 

1 1-77 

Vi+i = -V] —Vj + H~ 2axj, (59) 



12 



which has the sokition 



1 , l-V 



i-1 



i=0 



V 



(60) 



Substituting the expression for i/j, Eqs. ([5^ and (IMl) . 
imposing the periodicity condition yo = one finds that 



Vq 



2ari 



(l + 77)(l-77-) 



(61) 



k=0 



which can also be rewritten as 



yo = 



2ar] 



(l + ?/)(l-7f) 



E^'(^fc-X-i-fe)- (62) 



k=0 



We can calculate again the remaining jjj from the cyclic 
permutation property. With as given above, we define 
y[x] such that using Eq. (pT|) we have jjo = y[x]- Then 



or explicitly, 



2ar] 



(l + 77)(l-77^) 



(63) 
(64) 



mods X(s— 1 



+j — /c) mo 



d s} 



fc=0 



It can be shown that in terms of the hull- function f{x) 
of Aubry's solution Eq. (|6]), the above equation upon 
substitution of Eq. ([51]) reduces to 

ij, = / (j^ 2a) - 2a Int (j^) = .g (j^ 2a) , (65) 

with g{x + 2a) = g{x), from which Aubry's result^^ fol- 
lows under the identification 



(66) 



Vj = V] + 2a Int [j 



Once the lowest energy configuration has been found, 
the mode-locking intervals of /i over which a given av- 
erage spacing €/2a = r/s is a lowest energy configura- 
tion as well as other properties pertaining to this con- 
figuration such as the Peierls-Nabarro barrier are readily 
obtained^-. 

In other words, these are properties that can be ob- 
tained from the global minimizers alone. The global re- 
gion is terminated by shocks, and thus as long as the 
location of the shocks are not known, we do not know 
the extent of the global region. It is therefore impossi- 
ble to know which one-sided minimizers will flow towards 
which global minimizer. We thus turn next to the cal- 
culation of the flow pattern which will also allow us to 
understand further what happens at the boundaries of 
the mode-locking intervals. We restrict the analysis to 
the case s = 1, which already contains the relevant fea- 
tures of the general case. 



F. Steady-state Flow Pattern for s = 1 

The steady state flow pattern contains a single shock 
tree and we will consider only the case s — r ~ 1, for 
which xo = 1— ■ As we have pointed out before, for a 
single shock tree the global segment is always intersected 
by a shock insertion. As can be seen from the right panel 
of Fig. [TUl the particle insertion steps cause a bifurcation 
of the global intercept v into a left and right intercept. 
For the case xo = Ij the equation for hj that governs the 
evolution of the global intercept, Eq. becomes 



-J> + 2aj, 



(67) 



which implies that the global intercept v remains on the 
right half of the bifurcation, cf. Eq. (|4T|) . This is equiva- 
lent to saying that the corresponding shock tree is a right 
tree. From Eq. ([5^ we find 



-(2a -/i) 



(68) 



The left branch of the bifurcation must obey the re- 
cursion for v, Eq. (|B2[) . with — bj — 2a and we thus 
have 



7^iy'f + il-rj){b,~2a). 



(69) 



With reference to the right panel of Fig. [TUl consider 
the time t = Q~ at which the global intercept bifurcates. 
Denote the common ancestor at t = — r as v^l so that 
i^-i = '^0 + 2a. Following the secondary intercept for 
times jr with j > 0, moving into the cell coordinates 
and solving the recursion for we find that 



- (2a - n) 



1-r; 



2ary 



(70) 



For reasons that will be apparent soon, we label the 

~(k) 

intercepts of a shock tree at time t — a.s , with 
k — 0, . . . K so that the number of secondary intercepts 
and shocks is each k -I- 1 , as shown in the right panel of 
Fig. [TUj With this labeling ^g"^ denotes the newly in- 
serted shock at t = 0. From the time periodicity t of the 
flow pattern it then follows that 



and thus 



~(fc) 



~{k) = AO) 



(2a - fi) 



2ar] 



k+i 



(71) 



(72) 



Therefore in the unit cell coordinates, the time evolution 
of a secondary intercept when projected back onto the 
set of intercepts at t = is simply a shift v'^^ — v'^^^'' . 
This can be clearly seen in the right panel of Fig. [TOl 
where the open circles mark the origin of the unit cell. 

For a steady state flow pattern with k + 2 intercepts, 
the global intercept must map into itself, the secondary 
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FIG. 10: (Color online) Steady-state shock trajectories (left) and evolution of intercepts v (right) for ^/2a = 0.9 and Ao = 0.2, 
corresponding to i?/2a = 1/1. There are k + 2 = 4 shocks right after insertion. Left: The shock tree is of right type. The 
global shock is shown in dark blue, while secondary shocks are light blue. Right: Global intercepts are shown in red, while 
local intercepts are shown in green. At a shock insertion, a segment is bisected giving rise to a bifurcation of the v intercept 
associated with that segment. The resulting pair of intercepts has been linked to the parent intercept by dotted lines and the 
times t — TIT at which the bifurcation occurs has been offset to a slightly earlier time for clarity. See text for details on the 
labeling of shocks and intercepts. 



intercept k — k must disappear upon further evolution 
and a new intercept is created at shock insertion. Thus 
the number of secondary intercepts remains the same 
from one shock insertion to the next. Since intercepts can 
only disappear if their corresponding shocks merge with 
other shocks, a steady state of the shift pattern requires 
that the secondary shock associated with the segment 
k = K collides with the global shock. We have already 
found the location of the global intercept Pq. In fact, 
note that from Eq. ([7^ we have that i>Q°°^ — The 
of the local intercepts under time 



r,(fc) 







(fc+i) 

^0 



shift I/, 

evolutions then suggests to label the global intercepts 
and shock as i/q""' and iQ°°\ respectively. The labeling 
of intercepts and their corresponding shocks is shown in 

Fig.nni 

Thus given the steady state profile in the cell coordi- 
nate system, the ordering of the corresponding intercepts 
is 



~(oo) ~(k) ~(k-1) 



^0 



^(1) ^(0) 



with 



~(oo) 



(73) 



(74) 



The number of shocks constituting a shock tree, k + 2, at 
a given insertion time is directly related to the lifetime 



6tc of a shock from its insertion to its absorption by the 
global shock as k = Int ((Jtc/r). 

The locations of the secondary shocks can now be 
found as follows. Consider time t = and let 
■ ■ ■ 5 Co'''' denote the initial locations of the sec- 
ondary shocks with the labeling in correspondence with 
that of the associated secondary intercepts, as shown in 
the left panel of Fig. [TOl The subsequent positions at 
time t = jr are obtained from Eq. (|B4p . In terms of the 
:(fe) _ Jk) 



cell coordinates C] — Q — bj, we find 



1 1 

J I ^ J 



~(k-l) 



(75) 



where v\ ^' 



~(oo) 



2a which compensates for the 
wrapping around the cell boundary for fc = 0. The 
shift property of the secondary intercepts necessarily ap- 
plies to their associated secondary shocks as well, so 
that — S^l'^^'^. Furthermore, by the shift property 



s!^^ = Co^'" for all k,jk + j < k, so that without loss of 
generality we can restrict ourselves to j = 0. For a right 
tree, it is convenient to let the location of the newly in- 
serted shock in the cell coordinates be ^q"'' = +a and one 
finds 



so 



(k) 



arj 



0,1,2,. 



(76) 
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The global shock ^q""-* and k are still undetermined, 
since so far we have not dealt with the collisions that must 
necessarily occur. We have obtained all intercepts as well 
as the positions of the secondary shocks. The steady- 
state shift motion of the intercepts described above im- 
plies that secondary shocks should not collide with each 
other during their time evolution. We will now show 
this explicitly by determining the time required for two 
adjacent secondary shocks to collide. Due to the time pe- 
riodicity T it is sufficient to do the calculation at t = 0"*". 
The velocity of a secondary shock fp*^^ is given as 



^c^) = a; e^'^) 



(77) 



for fc = 0,1,2,. ..,k;, where due to wrapping around the 
unit cell boundary we have again i>Q 



2a. 



Letting tc'^'' be the time of collision between shocks k and 
fc + 1, and noting from Eq. ([55)1 that \*^t = (1 - r])/r], 



^-V Co 



(fc) 



?0 



V 



j{k) _ y{k+l) ■ 



(78) 



We find for A; < k, that t'^c^ /t = 1/(1 — 77) thus 
all secondary shocks will collide simultaneously, if at all. 
However, since ry € [0,1], tc > r, two secondary shocks 
cannot collide during the time evolution < t < r (ex- 
cept for the case 77 = 0, corresponding to an infinitely 
strong external potential which we will ignore). Since 
the flow pattern has time periodicity r, this moreover 
means that they can never collide. Thus the only colli- 
sion possible is between the global shock and its adjacent 
secondary shock. 

The global shock location can be determined by use of 
the area constraint. 



u{x, t) dx ~ 0, 



(79) 



which follows from the continuity of the internal energy 
iJi„t(a; -t- 2a, r) = Hint{x, t). We find 



r(oo) _ J_ 



2a — II 



V 



1 



V 



rj 



K+i 



(80) 



A few points are worth noting. While the locations 
of the secondary shocks in the cell coordinates, Eq. ([75)) 
are independent of fj,, the location of the global shock 
does depend on /i. Moreover, the factor r]~'^ in front of 
the term in rectangular brackets will diverge as k — >■ 00 
unless the expression in the brackets vanishes sufficiently 
fast, which for each 77 puts a constraint on /i as a function 
of K, establishing thereby the region of /i and 77 values for 
which a steady-state flow pattern with s = r = 1 can be 
obtained. 

We next look at the conditions under which the global 
shock ^Q°°'' can collide with its adjacent secondary shock 
^g'"' , as required by the flow properties of the steady state 



solution. Denoting by t 
quirement is 



the time of collision, the re- 



< 



< 1. 



(81) 



Next, ti"^ is found using Eqs. ^ and ^ with 



(k.+ I) 



~(oo) 



and 
equivalent to 



(n+l) _ r(oo) 



Eq. dH]) then turns out to be 



1 



[1 



77 



< 2a — fi < a 



1 



1 



V 



77 



2(k+1) 



(82) 

The disjoint open intervals defined above, together with 
their closure points, cover the range < 2a — /i < 
a(l — 77) /(I -|- 77), which is precisely the mode- locking 
region for the steady state fiows with s = r = l^i for 
a given rj. In Ref.'^^ this interval was calculated as the 
range of values of fj, for which the energy per particle 
in the lowest energy configuration with a given average 
spacing is minimum. In the present framework however, 
this interval arises from a restriction on the form of the 
fiow pattern at steady-state. 

Immediately to the left (right) of the points /i^ defined 



as 



2a - /Ik 



1 



1+7/ 



[1 



V 



(83) 



the steady-state profile contains one more (less) sec- 
ondary shock, while at /i^ the global shock and its adja- 
cent secondary shock collide at the time of the shock in- 
sertion. Within each open interval of Eq. ([82]) the steady 
state flow pattern has k + 2 shocks. We thus see from 
Eq. ([80]) that as 2a — /i approaches the phase boundary 
2a — /ioo — a(l — '7)/(l + ??), there is an accumulation 
of inflnitely many shocks at the global shock whose loca- 
tion ^0°°-* 0. Notice that this is also the location of the 
particles in the corresponding lowest energy conflgura- 
tion, confirming the result that at the phase boundaries 
the trajectory of the global shock coincides with that of 
the global minimizer^^. 

Finally, let us obtain from the inequality Eq. (1821) a 



bound on the location of the global shock 
result is 



00) 



a77 



K + l 



The 



(84) 



Note that the boundaries of the above intervals are the 
possible locations of secondary shocks, Eq. ([75]) . 



G. One-Sided Minimizers and Discommensurations 

The fiow pattern also reveals what happens if we pick 
a particular time t and look at the configurations of the 
semi infinite chain as the endpoint x moves from —a to 
a. The locations of the particles with respect to the unit 
cell can be read-off by noting where the corresponding 
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(1) (2) (3) (4) 




X 



FIG. 11: (Color online) Sample particle configurations of a semi-infinite chain for ^/2a = 0.8945 and Ao = 0.4, corresponding to 
l/2a = 1/1. Left: The light green line corresponds to the global minimizer, the lowest energy configuration. The darker green 
lines are a subset of one-sided minimizers. These correspond to lowest energy configurations of the semi-infinite chain with the 
position of the outermost particle fixed. Particle locations are shown as black circles. Global and secondary shocks are shown 
in dark and light blue, respectively. Right: The particle configurations relative to the external potential for the corresponding 
one sided minimizers on the left panel labeled as (1) - (4). Notice the presence of discommensurations in (3) and (4), where the 
unit cell contains an additional particle. The alternate coloring of the particles shows that when a discommensuration occurs, 
all the particles to the left of the discommensuration must have moved by a period 2a. The discommensurations occur precisely 
when the backwards flow of the characteristics changes from the right to the left half of the unit cell. The characteristics shown 
in red terminate at a newly inserted shock, they are known as pre-shocks. 



location on the one-sided minimizer lies with respect to 
the well minima (red circles) and the well boundaries (red 
triangles). Sample one-sided minimizers and the corre- 
sponding particle configurations are shown in Fig. Illl for 
the one-periodic case s = r — 1. 

As we move into the chain, t ^ —oo, the coordinate of 
the corresponding particle in the external frame has to 
decrease. For the one-sided minimizers, i.e. the flow of 
characteristics traced backwards in time, this means that 
whenever the particle location relative to the unit cell in- 
creases, it must necessarily have changed wells. However, 
a decrease in the coordinate, in particular a move from 
the right half (0, a) to the left half (— a,0) of the unit 
cell, implies that the particle remains in the same well 
and creates a discommensuration in the case s = r = 1. 

The case of flow patterns with multiple shock trees is 
similar. In general, by keeping track of the number of 
insertions after which a change of unit cell occurs, one 
can determine if the given potential well contains ad- 
ditional or missing particles, i.e. discommensurations. 
One-sided minimizers flowing in the (global) regions be- 
tween the shock trees asymptotically approach the low- 
est energy particle configurations without incurring any 
discommensurations. On the other hand, one-sided min- 



imizers starting in a region between the branches of a 
shock tree, upon entering the global region in between the 
trees, generate discommensurations. The boundary sep- 
arating these two type of regions is a pre-shock^^ , a char- 
acteristic that evolves into a shock, such as the dashed 
red lines in Fig. [TTl 

The life-time of a newly inserted shock, namely the 
number of insertions before it merges with the global 
shock, also corresponds to the maximum number of par- 
ticles counting from the end point of the chain within 
which a discommensuration can occur. As we have shown 
for the period-one case, the phase boundary of the corre- 
sponding domain in the fJ. — plane is marked by an 
infinite number of secondary shocks. This turns out 
to be true for the phase boundaries of all the domains 
£/2a = r/s— and implies that there are particle configu- 
rations with discommensurations arbitrarily deep inside 
the chain. 

We now turn to the calculation of the particle con- 
figurations associated with the minimizers for the case 
s = r = 1. Denote by 1^ the regions bounded by the 
secondary shocks, Ik = iio^Kio''^^) for fc = 1, 2, . . . , 
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which using Eq. (|76)) is given by 

/fe = (ar/^a?7'=-l). (85) 

Likewise denote by I^+i = (^o°°\ Co""*) the region 
bounded by the global shock and the left-most secondary 
shock. The one-r backwards flow maps Ik into Ik-i- Let 
the initial time be to = 0. The interval Ij- is terminated 
on its right end by le^'"^) so that the segment of u is 
u'"^~^\y) = '"^'^""^•'(y, 0~), given as 

ui'~'\y) = Xl{y-Dt'^)-Xoy, (86) 

which noting that X*_t — X*^t — Xqt, X*_t — 1 — rj, and 
A!i_T = (1 — ?7)/77 can be rewritten as 

ui'~'\y)r=il^r^)(^y-^y (87) 

The At — T backwards flow of the characteristics maps 
■ijk e h into Xk-i as 

Xk-i ^Vk- u-^^^\yk)T. (88) 

Expressing Xk-i in the unit cell coordinates, as yk-i = 
Xk-i+2a — ji and using Eq. ([72]) we obtain the backwards 
recursion for the minimizer 

yk-i=mk + M^-ri)ri^~\ (89) 
with yk & Ik- The solution is found as 

Vk-j = v'Vk + v'' iv^' - V') - "^aSjk- (90) 

1 + ?7 

This equation is valid for j — 0,1,2, k. The case 
j — k corresponds to the transition from the right to the 
left half of the unit cell. It can be shown that in order to 
bring the coordinate back into the unit cell an additional 
2a has to be subtracted, accounting for the last term in 
the above equation. 

Denote by Iq the interval between the left boundary of 
the unit cell and the global shock, 

Io^[-a,it'^). (91) 

As is apparent from Figs. [TT] and [SI for yo € Iq, it must 
be that yk € Iq for all fc < 0. We will now verify this 
explicitly. Notice that Iq is the interval belonging to the 
global segment of u. The corresponding intercept in the 

unit cell coordinates is given by £>^°°\ Eq. Working 
again in the unit cell coordinates, the backwards map for 
fc < turns out to be 

Vk-i = VVk (92) 
whose solution is given as 

y-k = v'^yo fc > 0. (93) 



It is clear that as fc — )■ oo, — monotonously, thus 
converging to the lowest energy configuration. Combin- 
ing Eqs. ([TO)) and we have thus explicitly shown 
that all one-sided minimizers converge to the global min- 
imizer y = 0. 

The coordinates t/j , J < fc of the corresponding config- 
uration in the fixed frame turn out to be given as 

r y, + (j - l)2a, 2<j<k 

yj = { vi, J- = o,i (94) 

[ yj -t- 2aj, j <0 

with ijk G Ik- The discommensuration is generated by 
the particles j — and 1, which are in the same cell. 

For a bi-infinite chain with particle fc fixed at yk 
such that the corresponding unit cell coordinate satisfies 
ijk & Ik with fc > 0, the corresponding configuration still 
contains only a single discommensuration. Note that the 
other semi-infinite half extending to the right is equiva- 
lent to a semi-infinite chain extending to the left with its 
end-point at —yk and the chain reflected around the axis 
y = 0. Due to the structure of the intervals Ik, it follows 
that if jjk € Ik with fc > then — y^ e Iq. Hence iijjk G Ik 
with fc > 0, the semi infinite chain extending to the right 
cannot contain any additional discommensuration. The 
bi-infinite chain contains thus a single discommensura- 
tion. 

There are also bi-infinite chain configurations with- 
out any discomniensurations. They are given by yo € 
i~io"\io°°^)- This is the region between the global 
shock and its image obtained upon reflection at y — 0. 
Note in particular that the lowest energy configuration 
generated from yo — belongs to this interval as well, as 
it should. 



V. DISCUSSION 

We have shown that the flow of characteristics asso- 
ciated with a forced inviscid Burgers equation is related 
to the lowest energy configurations of FK chains. The 
trajectories of these characteristics traced backwards in 
time are the minimizers: the one-sided minimizers gen- 
erate the lowest energy configuration of a semi-infinite 
chain for which the location of the outermost particle is 
fixed. They also converge to limiting trajectories, the 
global minimizers, which generate the lowest energy con- 
figurations of the bi-infinite chain. The flow of minimiz- 
ers is confined to channels bounded by the trajectories 
of shock discontinuities that emerge from Burgers evolu- 
tion. The shocks form tree-like structures and separate 
topologically distinct configurations of the semi-infinite 
chain that are marked by the presence or absence of dis- 
commensurations and their locations. The shocks and 
their evolution are a consequence of the weak solutions 
to Burgers equation and, as we have shown, follow from 
thermodynamical considerations. 

There are possible extensions of the approach pre- 
sented here. Flow patterns containing shocks imply that 
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the corresponding particle configurations are pinned by 
the external potential. In fact, the case of a piece- 
wise parabolic potential can be considered as the limit 
when the external potential is so strong that particles 
are mostly confined to the bottom of the potential wells, 
which can be approximated by parabolic segments. It is 
therefore natural to consider potentials that deviate from 
being piece-wise parabolic. As is clear from our results, 
the corresponding profiles will still consist of continu- 
ous segments terminated by shock discontinuities, but 
the segments will not be straight lines anymore and the 
flow pattern will be perturbed. It should be possible 
to carry out a perturbation calculation. Knowing where 
this might break down would shed further light on the 
relation between shapes of external potentials and the 
intricate phase diagrams for the structure of their lowest 
energy configurations. 

For steady-state flow patterns containing shocks, the 
corresponding Burgers profile will always contain a seg- 
ment that is bounded by shocks that will never merge 
and thus the segment will never disappear. In fact, the 
existence of such a global segment is guaranteed under 
rather general conditionsi^ii^. Since the global segment 
contains the global minimizers, we were able to obtain 
these by simply searching for the characteristic trajec- 
tories in this region having the appropriate periodicity. 
The flow in the global region also allows us to calculate 
the backwards flow of characteristics in the vicinity of 
the global minimizer to which they necessarily converge. 
However, as long as the location of the shocks marking 
the boundary of the global region are not known, it can- 
not be asserted whether these characteristics are genuine 
one-sided minimizers and thus correspond to lowest en- 
ergy configurations of the semi-infinite chain or not. To 
give an example, without the knowledge of the locations 
of shocks in Fig. 1111 the corresponding particle configu- 
rations on the right panel cannot be determined. 

Thus while it appears to be possible to calculate the 
global minimizers from limited local information of the 
flow, in order to calculate one-sided minimizers we re- 
quire the full flow pattern including shocks. This cor- 
responds to the limited extent of Aubry's theorem pre- 
scribing only the structure of the lowest energy config- 
urations associated with the global minimizers, but not 
those associated with the one-sided minimizers. For sys- 
tems with random external potentials, for which Aubry's 
theorem is not applicable and an exact analytical treat- 
ment might not be possible, one could still be able to de- 
termine the global minimizers even if only approximately. 
Such an approach is implicit in the work of Feigel'man^lj 
where a description similar to a periodically-forced Burg- 
ers equation was constructed for a charge-density wave 
system with random impurities with a focus on calculat- 
ing the effective impurity pinning strengths rather than 
the phase configurations. 

As we have shown, the description in terms of a peri- 
odically forced Burgers equation lends itself to including 
temperature, Eq. (PH)) . The evolution Eq. (|32p becomes 



now 

kT °° 
Ut + uux = — Uxx + ^^S{t- riT) V'{x + tit), 

n=0 

where u{x, t) is related to the free-energy e{x, t) of the 
semi-infinite chain as exix,t) = u{x,t). Note that for 
non-zero temperatures the viscous term smoothens out 
u(x,t). In the case of a piece- wise parabolic potential 
for which u{x, t) consists of linear segments, the primary 
effect of T will be a rounding of the discontinuities mark- 
ing its boundaries, while the interior of the segments will 
still remain approximately linear. Dimensional analysis 
shows that the length scale over which the shock discon- 
tinuity is smoothened out is of the order At ~ vfcTr. 
Thus one expects that regions with shock spacing of order 
At or less will coalesce. This can happen at the accu- 
mulation of shocks in the profile of u{x, t) near the phase 
boundary of a domain with a given ^, as well as when 
the trajectory of a minimizer flows close to a shock, such 
as the global minimizer in Fig. 1111 As we have seen, for 
both of these cases the distances are of the order 5 ~ ar]'^ . 
Thus when Ay and 5 are comparable we expect that the 
corresponding configurations will be susceptible to ther- 
mal fiuctuations, that can give rise to discommensura- 
tions. On the other hand, for those portions of the min- 
imizer that stay sufficiently far from shocks {5 ^ Ay) 
the segments of u remain still approximately linear, and 
they should therefore be less prone to thermal fluctua- 
tions. The discommensurations formed under thermal 
fluctuation were prescribed in=2, and are precisely of the 
form given in Eq. (j94p . This is what one would expect, if 
the temperature is sufficiently small so that the density 
of discommensurations is low. 
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Appendix A: Characteristic Flows and the Inviscid 
Burgers Equation 

In this appendix we briefly review the weak solutions 
of the inviscid Burgers equation. For a more detailed 
account see22r— . 

Eq. (jl6p is in the form of a hyperbolic conservation law 

^^^^ 

We are looking for a solution of 

Ut + uux = (A2) 
subject to the initial condition 

u{x,0) = uq{x). (A3) 
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Given Eq. (|A2p . we define its characteristics as the 
curves x{t) in the xt plane on which u{x, t) remains con- 
stant. These curves are straight hues given by the char- 
acteristic equation 



x{t) = + tuo{xa), 



(A4) 



with uq{xq) being the speed of the characteristic emerg- 
ing from the point xq. An imphcit solution is then found 



u{x,t) = uo(a;o), 



(A5) 



where for a given {x^t), xq is determined from the char- 
acteristic equation, Eq. (jA4p . 

Depending on the initial conditions, the characteris- 
tics can intersect, giving rise to multiple-valued points 
that are resolved by introducing discontinuities (shocks). 
Even with smooth initial data, discontinuities can de- 
velop in a finite time. Since solutions with discontinu- 
ities do not form a strict solution of the partial differen- 
tial equation, one denotes these as weak solutions which, 
instead of the local PDE, are required to obey a weaker 
form of the conservation law. 



dt 



dx x(x,t) 



du 
'dt 



du 
dx 



= 0, 



(A6) 



for any continuously differentiable function x(a;,i) with 
compact supporfe^i"— . This still does not uniquely de- 
termine the behavior of discontinuities. In general, this 
requires inspecting the microscopic evolution from which 
the continuum description arose. In the case of the mass- 
spring system of Section IIII A| the weak solutions follow 
from demanding that the internal energy Hintix,t) as a 
function of the end point of the spring is continuous. 

Given a discontinuous segment of u with the disconti- 
nuity at Xq , the speed of the characteristics immediately 
to the left and right of xq are given as ui = u{xq ,t) and 
Ur = u{xQ,t), respectively. There are two cases that one 
needs to distinguish: (i) ui > Ur, and (ii) ui < Ur- In the 
former case, we have a moving shock discontinuity, while 
in the latter case, we have a rarefaction wave. 

(i) ui > Ur'. Applying the integral form of the conser- 
vation law around the discontinuity, it can be shown that 
the shock moves with a speed 



V = ^{ui+Ur 



(A7) 



This is known as the Rankine-Hugoniot jump condition. 

(ii) ui < Ur'. In this case the characteristics imme- 
diately to the left and right of the discontinuity at xq 
diverge from each other. The weak solution in this case 
turns out to be given by 



u{x, t) = Ul + 



a; — xo 
t 



for (a;, t) such that 



< 



< Ur 



Ul, 



(AS) 



(A9) 



1. Shock Motion and Collisions 

For the FK model with piece-wise parabolic potential, 
u{x,t) is a series of straight line segments of identical 
slopes and discontinuities, as shown in Fig. H) Consider 
first the evolution of a single straight line with initial 
slope Aq > and x intercept z/, so that 



ua{x) = Xo{x - v). 



(AlO) 



Characteristic lines emerging from xq move towards the 
left (right) for xq < v [xq > ly) while the characteristic 
line emerging from xq — v remains stationary. From the 
characteristic equation we thus find that 



u{x, t) — \{t) {x — v) 



with 



Ao 



1 + Aot 



(All) 



(A12) 



If instead of a straight line we consider a line segment 
initially bounded by xi < Xr- the evolution of this seg- 
ment will again be given by Eq. (jAllI) with the restriction 
xi + uo{xi)t <x< Xr + uo{Xr)t and regardless of whether 
V lies inside or outside the interval bounded by xi and 

We consider next the motion of a single shock initially 
at ^0 such that the slopes of the segments immediately 
to its left and right are given by Aq > 0. Let and v' 
denote the locations where the segments to the left and 
right of the discontinuity intersect the x-axis. In order 
to have a shock discontinuity we also require that y < v' 
and uq{x) is given as 



uq{x) 



Xo{x ~ y) X < 
Xo{x -v') X > ^0 



(A13) 



The values of u, immediately to the left and right of the 
shock are 



Ul = Ao(<^o -v)>Ur = Ao(Co - v') 



(A14) 



and the initial speed of the shock is thus given by 
Eq. ISl 



= Ao Co 



v + v' 



(A15) 



From the definition of the shock speed Eq. (jA7|) . it 
is also clear that ui > v > Ur- meaning that as time 
goes on, characteristic trajectories in the left and right 
vicinity of the shock will collide with the moving shock. 
For any given time t, the characteristics that have not yet 
collided with the shock will evolve their associated line 
segments according to Eq. (jAlip . As we have seen above, 
this evolution is such that the interception points v and 
J/' remain stationary. The slope of these segments will 
be given by Eq. (|A12[) and denoting by £^{t) and v{t) the 
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position and velocity of the shock at time t, respectively, 
we find that 



v + v' 



(A16) 



Differentiation of v with respect to t gives w = 0, so that 
the shock moves at constant speed. Fig. [T^ shows the 
evolution of u{x,t) at three subsequent times tQ,ti and 
t2 along with the flow of characteristics and the trajectory 
of the shock. In terms of the characteristic flow a shock 
acts like an attractor, gradually absorbing characteristics 
along with the associated values of uq (x) that flow with 
them. 

Consider now two shocks moving towards each other. 
We will denote the shocks as £,i and ^r- The corre- 
sponding initial profile uq(x) will consist of three line 
segments, whose corresponding intercepts we will label 
as I'l < Vm < i^r- The initial profile is thus given as 



Ao(a; 

uq{x) = { \o[x 
\q{x 



vi) X <£,i 

^ra) £.1 <X <£,r 
Vr) X > 



while the corresponding speeds of the shocks are 



vi = Ao y^i 

= Ao £,r 



2 

Vrn + Vr 



and 



(A17) 

(A18) 
(A19) 



In order for the shocks to collide we must have that 



n - Vr 



> 0. 



(A20) 



Let tc denote the time of collision. As t approaches tc 
the segment between the two shocks narrows until it dis- 
appears at tc, leaving a single shock. The weak solution 
prescribes that for t > tc this shock will continue to move 
as a single shock with shock velocity v' . By applying the 
Rankine-Hugoniot condition Eq. (jA7p at the instant the 
two shocks have just merged into a single shock, one finds 
that v' is given by 



(l^m - l^l)vi + {fr - l^m)Vr = {f,- - l^l)v' . 



(A21) 



The above equation resembles the conservation of mo- 
mentum, and thus the two shocks behave like particles 
with (constant) "masses" {vm — vi) and {vr — Vm) that 
collide inelastically. 



Appendix B: Burgers Evolution of the FK Model 
with Piece-wise Parabolic Potentials 



The variables v^'^'' and ^'^'^^ along with the asymptotic 
value of the profile slope A'j_ completely determine u{x,t). 
During the evolution from jr < t < (j + l)r, the variables 

Uj''^ associated with segments k remain constant or dis- 
appear due to merger of shocks, while the non-colliding 
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FIG. 12: (Color online) (a) Time evolution of the profile 
u{x,t) containing a single shock. The intercepts of the left 
and right segments with the x axis are denoted as f and v' . 
The position of the shock discontinuities at subsequent times 
to,ti and t2 are labeled by ^o,Ci s-nd ^2. (b) Characteristic 
lines associated with the profile (in red) and world-line tra- 
jectory of the shock discontinuity (in blue). Note how with 
increasing time more and more characteristics merge with the 
shock. 



shocks evolve according to Eq. (|45|) . During the particle 
addition step the shock locations remain unchanged 
as a new shock is inserted at ^ but the i/^'^^ variables 
are mapped according to Eqs. (j40|) . (j4l]) and (|42l) . while 
the parameter b indicating the location of the zero in- 
tercept evolves according to Eq. (l46t . Together with the 
rules of how to handle colliding shocks, we thus have a 
discrete dynamical system for the variables v and ^ that 
underlies the evolution of u(a:;, t) under the forced Burgers 
equation. Denoting by subscripts j the times t = (jt)"*" 
right after shock insertion, the evolution equations for the 
segments that do not disappear during a shock collision 
become 
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The above equations assume that the segments k are 
not involved in the coUision of shocks. A segment k will 
disappear during the time interval [j'r, (j ' + l)r) , if 

Avf^>0 and AefVA«f^<r, (B6) 

where Az;f ^ ^ ^f-^^ - vf and A^f ^ ^ ^ - 

After the collision, ^'^^'^ will continue to move with a new 

velocity s that has been worked out in lA II 



1. Feeding-Order of Newly Inserted Shocks and the 
Evolution of the Global Intercept 

The global intercepts lie in the strips bounded by shock 
trees and each of these strips contains one global mini- 
mizer. Thus at any time t — tit there are s locations of 
the global minimizers which correspond to the s topolog- 
ically distinct positions of the particles in the lowest en- 
ergy configuration. Let us denote these locations by y^, 
with a = 0, 1, 2, . . . , s — 1 and G [—a, a)- Thus jja are 
the locations of the particles in the external frame pro- 
jected back into the unit-cell by translations of 2a. The 
labeling is such that yo is the equilibrium configuration 
of the particle closest to the left boundary, y — —a, of the 
unit cell, yi refers to the particle in the lowest energy con- 
figuration immediately to its right, §2 denotes its nearest 
next neighbor to the right etc. The labeling a is a num- 
bering of the particles according to their positional order 
in the lowest energy configuration. Observe that unless 
r = 1 the sequence is not monotonously increasing, 
since the period of the configuration will comprise r unit 
cells, whereas {jja} are the locations projected back into 
a single unit cell. 

Now focus on a single global minimizer. The location 
of this minimizer at a time t = nr must correspond to one 
of the {y}, say ya- Note that this location also marks the 
position of the particle at the end point of a semi-infinite 
chain. At the next insertion time t = (n+l)T the location 
of the global minimizer in the unit cell must necessarily 
be that of the next particle in the periodic configuration. 



say yp. With the labeling convention given above, we 
have /3 — (a + 1) mod s. The same is true for all other 
global minimizers. Thus from one insertion time to the 
next, the position of each of the global minimizers cycles 
through the ordered set {j/q}. 

On the other hand, at any given insertion time the lo- 
cations of the s global minimizers are distinct and they 
form the set {ya}. Thus we can also order the set 
of {ija} according to proximity in the unit cell [—a, a). 
Let us assume that the ordering in this way is given as 
{yacVai,- ■ ■ , where ao, ai, . . . , a^-i is some per- 

mutation of 0, 1, . . . , s — 1. It is not difficult to convince 
oneself that the differences (a^ — a^+i) mod s must 
be identical: Given a time t = nr, the location of the 
shock just inserted, by definition also marks the 

left boundary of the unit cell. Thus yo defined above 
as the global minimizer closest to the left boundary is 
also closest to the new shock from the right. At time 
t = (n + 1)t a new shock is inserted at £,'^"^1 = Q""" - fi, 
cf. Eq. (IBSP and thus there is a corresponding global 
minimizer immediately to its right corresponding to yo 
at this new time. Thus when progressing in time, the 
location yo must cycle through the s global minimizers 
which we had labeled as a = 0, 1, 2, . . . , s — 1, at some 
earlier time Iq — nr. The uniform shift by —fi of the 
location of the new shock to be inserted implies that this 
cycling of yo through the minimizers must also be a shift 
of the form a — > (a — A) mod s, where A < s and 
A and s are co-prime. In fact, r = A, so that this can 
be regarded as a definition of r. Thus for a steady state 
fiow pattern corresponding to ^/2a — r/s, s determines 
the periodicity in time sr, while r controls the "feeding 
order" of the shock trees. 

Observe now that the feeding order of the shock trees 
also determines whether the shock associated with the 
right boundary of a global segment is to the immediate 
left or right of the newly inserted shock in the co-moving 
coordinates: Recall that (i) to each inserted shock there 
corresponds a shock tree into which this shock will even- 
tually fiow, and (ii) that for any any two neighbor- 
ing global minimizers are separated by a shock tree (and 
hence a global shock). The sequence of being to the left 
or right of the newly inserted shock must therefore also 
follow the feeding order. 

We thus find from Eqs. dHT]) and dE]) that for the 
global intercepts v^''^ on a shock tree 
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with each value oi 6 — 0, 1, 2, 
with one of the s shock trees. 



1 being associated 



^ Ya. Frenkel and T. Kontorova, Phys. Z. Sowtetumon 13, ^ M. Weiss and F. -J. Elmer, Phys. Rev. B 53, 7539 (1996). 
137 (1938). 



21 



^ D. Cule and T. Hwa, Phys. Rev. Lett. 77, 278 (1996). 

* N. I. Gershenzon, V. G. Bykov and G. Bambakidis, Phys. 

Rev. E 79, 056601 (2009). 
^ M. Peyrard, Nonlinearity 17, Rl 2004. 
® The Prenkel-Kontorova Model - Concepts, Methods and 

Applications, O. M. Braun and Y. S. Kivshar, Springer, 

Berlin Heidelberg, 2004. 
^ S. Aubry in Solitons and Condensed Matter Physics, cd. 

A. R. Bishop T. Schneider, Solid State Sciences 8 264, 

Springer, Berlin, 1978. 
® S. J. Shenker and L. P. Kadanoff, J. Stat. Phys. 27, 631 

(1982) . 

® S. N. Coppersmith and D. S. Fisher, Phys. Rev. B 28, 2566 

(1983) . 

^° M. Peyrard and S. Aubry, J. Phys. C: Solid State Phys. 
16, 79 (1983). 

" R. S. MacKay, Physica D 7, 283 (1983); R. S. MacKay, 
Physica D 50, 71 (1991). 
S. Aubry, Physica D 7, 240 (1983). 
J. N. Mather, Topology 21, 457 (1983). 
H. R. Jausslin, M. V. Kreiss and J. Moser, Proc. Symp. 
Pure Math. 65, 133 (1999). 

W. E, K. Khanin, A. Mazel and Ya. Sinai, Ann. Math. 
151, 877 (2000). 

W. E, Comm. Pure Appl. Math 52, 811 (1999). 
I'' A. N. Sobolevskii, Mat. Shomik 190, 1487 (1999). 
^® J. Bcc and K. Khanin, Phys. Rep. 447, 1 (2007). 

R. B. Griffiths, H. J. SchcUnhuber and H. Urbschat, Phys. 

Rev. B 56, 8623 (1997). 
^° S. -C. Lee and W. -J. Tzeng, Phys. Rev. B 66, 184108 

(2002). 

^1 S. Aubry, J. Phys. C: Solid State Phys. 16, 2497 (1983). 
Classical dynamics: a contemporary approach, J. V. Jose 
and E. J. Saletan, Cambridge University Press, Cambridge, 
1998. 



2^ S. Aubry S and P. Y. Le Daeron, Physica D 8, 381 (1983). 

J. M. Greene, J. Math. Phys. 20, 1183 (1979). 

J. D. Cole ,Commun. Pure Appl. Math. 3, 201 (1950). 
^® E. Hopf, Quart. Appl. Math. 9, 225 (1951). 

Numerical Methods and Conservation Laws, R. J. LeVeque, 

Birkhauser Verlag, Basel, 1992. 

Linear and Nonlinear Waves, G. B. Whitham, Wiley, New 
York, 1974. 

Partial Differential Equations, L. C. Evans, American 
Mathematical Society, Providence, 1998. 
^° R. B. Griffiths and W. Chou, Phys. Rev. Lett 56, 1929 
(1986). 

W. Chou and R. B. Griffiths, Phys. Rev. B34, 6219 (1986). 
In fact as we will see shortly, for sufficiently large times 
t, the flow reaches a steady-state limit in which the initial 
condition docs not matter anymore. 
T. Tatsumi and S. Kida, J. Fluid Mech. 55, 659 (1972). 
Only in the special case where at a given insertion time t 
the shock tree contains only one shock, namely the global 
shock, the corresponding global intercept will also be the 
one associated with the global shock. 
M. Mungan in preparation. 

^® The case for r = and s = 1 so that xo = is analogous, 
but unphysical from a classical physics point of view, since 
the average spacing of particles in the lowest energy con- 
figuration vanishes. This results in a point-like condensate 
at the minimum of the potential well. 
M. V. Feigcl'man, Sov. Phys. JETP 52, 555 (1980). 
Controlled Markov processes and viscosity solutions, 
W. Fleming, H. M. Soner, Springer, New York, 1993. 

^® D. Gomes, R. Iturriaga, K. Khanin and P. Padilla, Moscow 
Math. J. 5, 613 (2005). 

F. Vallet, R. Schilling and S. Aubry, J. Phys. C: Solid State 
Phys. 21, 67 (1988). 



